An Erupting Classical Nova in a Globular Cluster of M87 
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ABSTRACT 

Only one certain classical nova eruption has ever been detected inside a glob- 
ular cluster - nova 1860 A.D. (T Sco) in M80. During a survey of M87 we 
have detected an erupting star coincident (to within 0.08 pixels) with a globu- 
lar cluster of that giant elliptical galaxy. We are able to discount variables in 
the foreground or background of M87. The light curve and color of the erupt- 
ing star match those expected for a nova at the distance of M87. The chance 
superposition of an M87 field nova on the globular cluster is very unlikely but 
cannot be completely ruled out. Our detection hints at a globular cluster nova 
frequency / ~ .004 novae/cluster/year, much higher than previous observations 
have suggested. 

Subject headings: clusters: globular; cataclysmics: novae 

1. Introduction 

The realization that X-ray sources are a thousand times more common in globular 
clusters than in the rest of the Galaxy (Clark (1975)) first spurred theorists to consider the 
effects of very close stellar encounters. Tidal captures may occur during very near-misses of 
degenerate stars (Fabian, Pringle & Rees (1975)) and low mass main sequence stars. Neutron 
stars with captured companions give rise to the observed cluster X-ray sources, but white 
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dwarfs are much more common than neutron stars. Thus dozens of tidal-capture white dwarf 
- red dwarf binaries (cataclysmic variables, including classical novae = CN) are also expected 
in dense globular cluster cores (Rappaport &di Stefano (1994)). Nearly three decades of 
predictions (e.g. Fabian, Pringle & Rees (1975), Krolik (1984)) and surveys (e.g. Grindlay 
et al (2001), Knigge et al (2002)) have finally begun to reveal large numbers of cataclysmic 
binaries (CVs) concentrated in the cores of Galactic globular clusters. Contrasting tidal 
capture-generated novae with binary evolution-generated novae would be very interesting. 
For example, do the two populations of novae outburst with similar frequencies? It is clearly 
necessary to first find at least a few tidal capture novae to begin such a study. 

Over the past century, of order a thousand erupting CN have been detected in the 
Galaxy and M31. However, only two candidate CN eruptions have ever been observed in 
Galactic globular clusters (GC) (Luther (1860) and Hogg & Wehlau (1964)). Of these only 
nova 1860 (T Sco in M 80) seems beyond any doubt to be a cluster member. This is based 
on the remarkable and independant discoveries and light curve documentation by Auwers 
(1862) and by Pogson (1860), supplemented by Sawyer (1938) and Wehlau, Sawyer-Hogg 
& Butterworth (1990). Shara & Drissen (1995) have recovered the quiescent nova in HST 
images of the core of M 80. The second candidate (nova 1938 in M 14) was a genuine 
variable that had the right brightness to be a CN, but details of its light curve are far too 
scant to determine its true nature (Hogg & Wehlau (1964)). Recovery attempts have been 
unsuccessful (Shara et al (1986); Margon et al (1991)). 

If classical novae were as overabundant in globular clusters (relative to the field) as are 
X-ray sources we should have detected about 100 CN in the GC of the Galaxy and M31 over 
the past century. While there might be one or more CN depletion mechanisms in GC, the 
present study suggests instead that searches for CN in GC are very incomplete. 

In the course of an HST imaging survey of M87 effectively spanning ~ 90 days we have 
recently detected an erupting CN coincident with a GC of that galaxy. If the nova rate in 
M87 GC is similar to that in the Galaxy we would expect to observe a nova only once in 
~ 20 years of HST observations (see section 4). The rarity of detected Galactic and M31 
GC CN, the ease with which the M87 nova was detected, and the possible imphcation for 
globular cluster nova rates prompt us to report it here; full details and further searches will 
be described elsewhere. 
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2. Observations Astrometry 

M87 was imaged by HST on the 30 successive days from 28 May 2001 through 25 June 
2001 inclusive. A log of the observations, taken through the F814W and F606W filters, is 
given in Table 1. We have also examined all the HST archival images of M87 and its globular 
clusters, so that our effective survey time is about 90 days (see below). Most other HST 
images of M87 and its globular clusters are single-epoch. Erupting classical novae brighten 
from M ~ +4 to M ~ — 8 in about one day. We would have thus detected the nova reported 
here in a similar, average luminosity M87 GC at maximum or within the ~ 5 days following 
maximum light. Thus the effective survey time for the 10 other epochs is ~ 6 days each, and 
for the entire survey it is ~ 30 + 60 = 90 days. The M87 GC in which we detected a variable 
is marked in the finder chart shown in Figure 1. The GC is located at q;(2ooo) =12:30:45.4, 
5(2000) =12:23: 16.1 as determined using the world coordinate system in the WFPC2 header. 
The GC is located 13.51 arcseconds south and 60.95 arcseconds west of the M87 core. While 
the GC is quite small and faint, we have determined that it is a resolved object and not a 
star. The globular cluster appears at / = 23.06, corresponding to Mj — —8.06 ... near the 
mean luminosity of globular clusters. The cluster displays V — I — +1.0 , also typical of 
globular clusters. 

In Figure 2 we show a series of five images centered on the position noted above. We 
dehberately added, at the left of Figure 2, an image taken six years prior to the others. This 

(and other epoch images) clearly shows the presence of the M87 GC without the variable. The 
circled GC, coincident with the variable, is close to, and significantly fainter than another 
GC above and to the right in Figure 2. The middle image of Figure 2 (8 June 2001) shows 
the variable at maximum light: it has effectively brightened its host GC to match the nearby 
GC. Towards the end of the 30 day observing campaign the variable has faded from view 
and the circled GC is again faint compared to its neighbor. 

The light curve of the GC plus variable is shown in Figure 3. The hint of a brightness 
increase on day 11 is strongly confirmed by the brightness peak seen in both F814W and 
F606W on day 12. The GC remains obviously brighter than its quiesent state for ~ 6 
days after maximum hght in the F606W filter and ~ 3 days in the F814W images. The 
photometry of Figure 3 is hsted in Table 1. 

3. Is the Variable a Nova? 

Have we discovered an erupting nova in an M87 GC? We show below that the light curve 
is fully consistent with a CN in M87. However, to answer the question in the affirmative we 
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must first rule out other possibilities. Most likely amongst these are: a chance superposition 
of a background supernova, a gamma ray burst, an M87 nova, a microlensing event or a Milky 
Way variable along the line of sight to the GC. The supernova hypothesis can immediately 
be discarded because of the rapid rise and decline times of the light curve. 

The detailed shape of the decline from maximum light on day 12 through day 21 suggests 
a power law decline L a t~^'^ . There are observed optical afterglows of gamma ray bursts 
(GRB) with similar dechne laws (e.g. GRB000301C Jensen et al (2001) and GRB980703 
Holland et al (2001)). GRB afterglows typically achieve R or I ~ 20 one day after outburst; 
the variable in M87 is 3 magnitudes fainter. One might thus argue that the outbursting 
object is at higher redshift than most GRB afterglows observed so far. We regard this 
explanation as very unlikely both because of the rarity of GRB (~ 1/day over the entire 
sky) and because there is some indication of a precursor rise in the F606W image on day 11. 

The constant brightness of the GC in images previous and subsequent to the eruption, 
the rapid brightness rise and several day decline, and the brightness and moderately blue 
color of the variable rule out most types of Galactic variables. RR Lyraes, eclipsing binaries 
and flare stars don't match the observed brightness profile and/or blue color. If the variable 
is a dwarf nova along the hue of sight to M87 it is about 100 kpc distant from the Galaxy... an 
intergalactic tramp. We regard this possibility as extremely unlikely. The asymmetry of the 
light curve likely rules out a microlensing event. 

In F606W the brightening of the GC is ~ 1.2 magnitudes. Taking into account the 
brightness of the GC the peak mag of the variable in F606W is 23.57. At the distance of 
M87 (m — M) ~ 31.12 (Kundu et al (1999)) this corresponds to an absolute magnitude of 
—7.55. The slope of the variable's decline yields decline times of ^2 ~ 18 days and ^3 ~ 27 
days. Using the absolute mag - decline time relations for novae from Downes & Duerbeck 
(2000) yields My ~ —8.35 (ts relation) and My ~ —8.12 (t2 relation), in good agreement 
with the observed value. The moderately blue color of the variable (at peak brightness it 
displays V - I = 0.61) further supports the suggestion of a classical nova. 

The most difficult possibility to rule out is the chance superposition of an M87 field 
nova on the GC. We have been quite successful in detecting erupting novae in M87 (Shara 
& Zurek (2002a) and Shara & Zurek (2002b)). About 3 to 5 CN in eruption (~ 1 or 2 on 
each WF chip) are detectable with HST in F814W band images every time one looks. Novae 
are usually much brighter in ultraviolet or blue passbands so the number in eruption in each 
epoch may be considerably higher. This allows us to crudely estimate the probabihty of 
a chance GC-CN superposition. Each WF chip has 562, 500 useful pixels, and about 0.5% 
are covered by M87 GCs. Randomly distributing 1 nova on a WF chip leads to a chance 
cluster-nova miss rate of 0.995 and thus a chance coincidence rate of 1 — 0.995 = 0.005. 
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The near perfect positional coincidence of the GC and the nova candidate (see Figure 4; 
±0.08 pixels, corresponding to 0.6 pc at the distance of M87) make the chance coincidence 
rate ~ 100 times better ... i.e. about 0.0005%. We understand the dangers of a posteriori 
statistics, and recognize that we may well have detected a fluke superposition of a field CN 
on an M87 GC. In the rest of this letter, though, we will assume that the CN is a genuine 
globular cluster member. 

4. Nova Rates 

The detection of one certain erupting CN in the ~ 150 GC of the Milky Way over the 
past ~ 140 years places a rough lower limit on the nova-in-cluster eruption frequency f: 
/ > .00005 CN/cluster/year. Weekly observations are required to find the fastest novae, but 
most Galactic globular clusters have not been observed this frequently (by either professional 
or amateur astronomers). Erupting CN rival their host globular clusters in luminosity and 
should draw attention to themselves, but only to amateur and professional astronomers 
realizing that a very bright star in a globular cluster is unusual. The incompleteness factor 
in f is difficult to estimate, but it must be at least a factor of 2 (due to seasonal and weather 
effects) and is probably much more. 

Motivated by theory and models suggesting the existence of many cataclysmic variables 
(CVs) in GCs, CiarduUo, Tamblyn & PhiUips (1990) and Tomaney, Crotts & Shaffer (1992) 
surveyed large numbers of M31 globulars for erupting CN. None were detected in an Ho; 
survey over an effective survey time of one year for the entire known M31 GC system. This 
places a useful upper limit on f: / < .005 CN/cluster/year. (We note in passing that while 
the range of observed globular cluster luminosities and masses is substantial, the overall 
Galaxy and M31 globular cluster systems are remarkably similar (Racine &Shara (1979)). 

The survey we are reporting here had an effective time coverage of 90 days for the 1057 
globular clusters of M87 reported by Kundu et al (1999). This one detected erupting nova 
in our survey yields a GC nova frequency f~ .004 CN/cluster/year, just slightly smaller than 
the upper limit reported by Tomaney, Crotts & Shatter (1992). Unless we were very lucky 
to have detected a GC CN (or the CN is not a GC member) this suggests that classical nova 
eruptions in globular clusters are up to 100 times more common than current detections in 
the Milky Way suggest. A search in M31 extended in time by a factor of about 3-5 over that 
of ?)tcs92) should reveal two or more GC novae, and firm up or refute the tentative (and 
surprisingly high) rate we report in this letter. While a twice- weekly search of every globular 
cluster in the Milky Way is a daunting task, it is not beyond the capabilities of modest 
robotic telescopes. A rate of .004 novae/yr/ cluster suggests a remarkable overall nova rate 
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of 1 event every 2 years in the globulars of the Milky Way... well worth searching for. 

We gratefully acknowledge support under grant GO 8592 from the Space Telescope 
Science Institute. This research is based on observations made with the NASA/ESA Hubble 
Space Telescope obtained at the Space Telescope Science Institute. STScI is operated by 
the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 
5-26555. 
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Fig. 1. — Finder chart of the erupting classical nova in M87 which is coincident with a 
globular cluster. The square in the lower left quadrant of the HST WFPC2 image is shown 
magnified in the overlay to its right, in which the GC host is indicated with an arrow. See 
figure 2 for a closeup. 
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Fig. 2. — Closeup of the erupting classical nova coincident with the globular cluster of M87. 
Left image is from a F555W WFPC2 frame taken 1995; the four rightmost images, taken 
in 2001, show the rise and decline of the nova (circled with its host globular). A second, 
brighter field globular cluster of M87 is above and to the right of the circled nova and GC. 
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Fig. 3. — F606W and F814W light curves of the M87 nova and globular cluster combined. 
Magnitudes are in the ST system. 
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Fig. 4. — Comparison of the globular cluster (center of each panel) without the nova erupt- 
ing (left) and with the nova erupting (right). The distribution of light during eruption is 
consistent with the nova being superposed directly on the GC. These subimages are taken 
from the interlaced F814W frames; the pixel scale is twice as fine as the F606W images 
shown in Figure 2. The left panel is the average of all F814W images that do not have the 
nova in outburst. The right panel is the F814W image from day 12. 
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Table 1. Observations and Photometry (30 days beginning 2001/05/28 and ending 

2001/06/25) 
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